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ABSTRACT 

We show that the recently observed elemental abundance pattern of the carbon-rich metal-poor 
Damped Lyman a (DLA) system is in excellent agreement with the nucleosynthesis yields of faint 
core-collapse supernovae of primordial stars. The observed abundance pattern is not consistent with 
the nucleosynthesis yields of pair-instability supernovae. The DLA abundance pattern is very similar 
to that of carbon-rich extremely metal-poor (EMP) stars, and the contributions from low-mass stars 
and/or binary effects should be very small in DLAs. This suggests that chemical enrichment by the 
first stars in the first galaxies is driven by core-collapse supernovae from ~ 20 — SOM© stars, and also 
supports the supernova scenario as the enrichment source of EMP stars in the Milky Way Galaxy. 
Subject headings: galaxies: abundances — galaxies: evolution — quasars: absorption lines — stars: 
abundances — stars: Population HI — stars: supernovae 



1. INTRODUCTION 

At the end of the dark age of the Universe, the cos- 
mic dawn was heralded by the birth of the first stars 
and galaxies. The nature of these first objects is still 
far from being well understood. From the theory of 
star formation from primordial gas, the first stars are 
believed to be very massive, with masses of the or- 
der of IOOM0, given the limited cooli ng of molecular 
hydrogen (e.g., lIBromm fc Larsonl[200l . This depends 
on fragmentation in a cosmological minihalo, ionization 
prior to the onset of gravitational collap se, and the 
accre t ion rate from the cl oud envelope CTan fc McKed 
120041: lOhkubo et all I2009D . and it seems possible to 
form lower-mass stars (~ 10 — AO M(t,) from primordial 
gas in recent num erical simulations (|Yoshida et al.l[2b08l : 
iStacv et"al] l2009f ). The primordial stars with initial 
masses of ~ 140 — 270Mf7^ explode as pair instability su- 
perno vae (PISNe) (IB arkat e t al.llT967HUmeda fc Nomotol 
120021: iHeger fc Woo slcv 200^ In terms of chemical 
abundances, however, no observational signature for the 
existence of PISNe has been detected. 

Chemical enrichment from the first generation of stars 
has mainly been studied with extremely metal-poor 
(EMP) stars (e.g., .Beers & Christlicb 2Q05|) . 10 - 25% 
of EMP stars with [Fe/H] < - 2 IjAokil I2010D show 

carbon enhancement relative to iron ([C/Fe] >^ 1); 

this is also the case for the three hyper/ultra metal- 
poor stars known, with [Fe/H] < —4.5. Such stars 
are collectively known as carbon-enhanced metal-poor 
(CEMP) stars, and are further subdivided as CEMP- 
s and CEMP-no stars depending on whether they ex- 
hibit or not an enhancement in the abundances of slow 
neutron capture elements such as Barium. The abun- 
dance patterns of EMP stars are explained with one 
of, or the combination of, four enrichment sources : 
(1) core-collap se supe rnovae (lUmeda fc: Nomotol 120031 : 
llwamoto et al.1.2005 : Tominaga et al.ll2007t ). (2) rotating 



massive stars (iMevnet et al.l[2006[ ). (3) asymptotic giant 
blanch (AGB) stars in binary systems (Suda ct al. 2004), 
or (4) interstellar accretion ((Yoshiil llQSlt llbeni I1983D . 
These enrichment sources have distinct signatures, so 
that it should be possible to distinguish between them 
by comparing the observed elemental abundances with 
theoretical calculations of nucleosynthetic yields. 

The high [C/Fe] can be explained with the mass trans- 
fer from AGB stars in binary systems, the mass loss 
from rotating massive stars, and a single core-collapse 
supernova forming a black hole. The Ba enhancement 
can be explained by AGB stars and p ossibly by rotat- 
ing massive stars (jPignatari et al.l[2008f ). However, iron- 
peak elements have to come from supernovae. The ef- 
fect of interstellar accretion has not been studied with 
hydrodynamical simula tions, but is obser vationally esti- 
mated to be negligible (jFrebel et al.ll2009| ). The binarity 
of CEMP stars has been studied by monitoring radial 
velocity variations. A signature of binar ity is seen sta- 
tistically for possibly all CEMP-s stars ( Lucatello et al.l 
|2005), but the bin ary fraction of CEMP -no stars seems 
to be much lower (|Aokill2010D . In particular, the bina- 
rity is not seen for the t hree stars with [Fe/H] < —4.5. 
lUmeda fc Nomotol (poM ) were the first to show that the 
observed abundance pattern from C to Zn can be well re- 
produced with the enrichment from a single core-collapse 
supernova that leaves behind a relatively massive black 
hole, under the assumption of inhom o geneous chemi- 
cal en richment (jAudouse fc Silk|[l995[ ). llwamoto et al.l 
(|2005t ) showed that the N abundance can be as large as 
observed with enhanced mixing between H and He layers 
during the hydrostatic stellar evolution. 

The observations of very metal-poor Damped Lyman 
a (DLA) systems have opened a new window to study 
the chemical enrichment of the Universe by the first 
generations of stars. DLAs are quasar absorbers de- 
fined by their high column density of neutral hydrogen, 
log A^(Hi)/cm~^ > 20.3. They appear to sample a range 
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of galaxy types, from the extended Hi disks of galaxies, 
to smaller subgalactic size hal oes, as well as sm aller Hi 
clouds within larger galaxies () Wolfe et alj 120051) . Large 
scale surveys, such as the Sloan Digital Sky Survey, have 
increased ten-fold the number of kn own PL As, which 
now number in excess of ~ 1000 (jNoterdaeme et all 
120091: iProchaska fc Wolfell2009l ). Follow-up high resolu- 
tion spectroscopy of the most metal-poor DLAs is of par- 
ticular interest, since the gas they trace m ay have been 
enriched by very few generations of stars (jPettini et alj 
[2008; Pcnprase et al. 2010). Moreover, measuring ele- 
mental abundances in DLAs is straightforward; the only 
potential complications are line saturation and dust de- 
pletion and both effects are of much reduced importance 
for metallicities Z < 1/100.^0. In addition, the chem- 
ical evolution of such systems is rather simple, whereas 
in more chemically evolved systems, there are uncertain- 
ties in the star formation history, gas inflow and outflow, 
large contributions from AGB stars and Type la super- 
novae, with the result that the signatures of the first stars 
can be easily washed out. 

Thus, if the most metal-poor DLAs hold the key to 
unravelling the chemical enrichment from the first stars, 
it is of great interest to compare their observed abun- 
dance pat terns with th e nucleo synthesis yields of metal- 
free stars. iCooke et al.l ()2010b[ ) recently reported such a 
DLA with [Fe/H] ~ —3, which exhibits a strong carbon 
enhancement relative to all other available elements, in- 
cluding [C/Fe] ~ -1-1.53. This reminds us of the CEMP 
stars in the solar neighborhood. In this letter, we com- 
pare the elemental abundance pattern of the C-rich DLA 
with the nucleosynthesis yields of both core-collapse and 
pair-instability supernovae (§2). In §3 we give a more 
general discussion of the chemical enrichment of the Uni- 
verse by the first generation of stars. We summarize our 
main conclusions in §4. 

2. ABUNDANCE PROFILING 

In stars with initial masses of > 10A/q, the iron core un- 
dergoes gravitational collapse at the end of the star's life. 
If sufficient amount of the released gravitational energy 
is transported to the outgoing shock wave, a supernova 
explosion occurs. The mechanism of core-collapse explo- 
sion and formation of a black hole remnant is still uncer- 
tain, although a few groups have presented feasible cal- 
culations of e xploding su pernovae (M arek fc Janka 20dl 
e.g.,) Thu s lUmeda fc Nomoto (2002I). Ilwamoto et all 
(|2005D . and iTominaga et al.. (,2007. ) have calculated nu- 
cleosynthesis yields by promptly exploding the pro- 
genitor model without calculating further collapse and 
bounce, and by applying mixing-fallback. During the 
supernova explosion, the elements synthesized in differ- 
ent stellar layers mix and a large fraction of this mixed 
material falls back onto the remnant black hole. Physi- 
cally, such a mixing-fallback process can be caused not 
only by the Ray leigh- Tay lor instability (jHachisu et al.l 
119901 : iJoggcrst et al. 2009j but also by a jet-induced ex- 
plosion (Tominaga 2009). In the present study, we cal- 
culate nucleosynthesis yields with two new models of 
faint core-collapse supernovae applying different mixing- 
fallback parameters. 

Figure [1] shows their elemental abundance ratios from 
C to Zn relative to Fe. We adopt the progenitor star 
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Fig. 1. — The elemental abundance pattern of the metal-poor C- 
rich DLA (filled circles) and peculiar DLA (open circles). The solid 
and short-dashed lines show the nucleosynthesis yields of faint core- 
collapse supernovae from 25Mq stars with mixing-fallback. The 
dotted line is for pair-instability supernovae from 170Mq stars. 

model from Ilwamoto et al.l (|2005l ) for an initial mass of 
M = 25Mq and metallicity Z = 0. We implement two 
cases with different explosion energies: 1 x 10^^ erg (su- 
pernova (SN), solid line) and 20 x 10^"'^ erg (hypernova 
(HN), short-dashed line). An efficient mixing- fallback 
is adopted in both models: the mixing region of ejecta 
is assumed to be Min — Mout = 1-76 — G.IAMq and 
2.20 - 6.30Mq with the ejection fraction of / = 0.008 
and 0.004, respectively. The resultant black hole masses 
are Q.IMq and 6.3Mq, respectively. The ejected Fe 
masses are as small as M(^^Ni) = M(Fe) = 0.0018^^0 
and O.OO14Af0, respectively, and thus correspond to faint 
supernovae such as SN1997D. For comparison, Af(Fe) ~ 
O.OTMq for normal supernovae and ^O.O7M0 for normal 

hypernovae in iKobavashi et al.l (hereafter K06, 120061 ). 
The parameter dependencies of elemental abundance ra- 
tios can be summarized as follows. 

• The abundances of odd-Z elements are lower than 
those of the a elements. This odd-even effect 
with the order of ~ 0.5 dex is realized for pri- 
mordial supernovae in general, and is consistent 
with observations of EM P stars in the Galactic halo 
(e.g., ICavrelet al.]|200l . The yields of odd-Z el- 
ements are enhanced by the surplus of neutrons 
in ^^Ne, and ^^Ne is transformed from ^''N during 
He-burning. The absence of seed CNO elements 
in the progenitor stars results in smaller amounts 
of odd-Z elements. This metallicity effect gives 
excellent agreement with the decreasing trends of 
[(Na,Al,Cu)/Fe] toward lower metallicity in the so- 
lar neighborhood (K06). The metallicity effect is 
stronger for Na, Al, Sc, and Cu than for P, CI, K, 
V, Mn, and Co. 

• C and N are synthesized in the outermost region 
of the ejecta. Thus, the [C/Fe] ratio is higher for 
smaller ejection of Fe. [N/Fe] depends on the initial 
metallicity and mixing during the stellar evolution. 

• a elements such as O, Mg, Si, and S are synthesized 
in the outer regions of the ejecta. [0/Fe] depends 
both on the explosion energy and mixing-fallback, 
and thus it is not possible to put a constraint on 
the explosion energy only from [0/Fe]. Since Mg 
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is more sensitive to the explosion energy, [Mg/Fe] 
may give a clue to solve this degeneracy. Si and S 
are synthesized in more inner regions than O, and 
significant amounts of Si and S fallback. 

• Zn, Co, and V are synthesized in the deep complete 
Si-burning region together with ~ 80 — 90% of Fe. 
[(Co,Zn)/Fe] > are realized in high energy ex- 
plosions with mixing-fallback (jUmeda fc Nomotol 
l2002f): ne utrino interactions could enhance [Zn/Fe] 
(|Pruet et al. 2005; Frohhch c^aD I2OO60 but not 
[Co/Fe] (jlzutani fc Umedall201C 



• Cr and Mn are mainly synthesized in the incom- 
plete Si-burning region. [Cr/Fe] does not vary very 
much for these models; Cr abundance does not de- 
pend very much on parameters such as the electron 
fraction and entropy, and [Cr/Fe] cannot be smaller 
with smaller Min or larger Mout- Larger Min can 
give larger [Cr/Fe] and smaller [(Co,Zn)/Fe], but 
such abundance pattern is not consistent with the 
EMP observations. Therefore, in galactic chem- 
ical evolution models, [Cr/Fe] is almost constant 
over a wide range of [Fe/H]. This has been a prob- 
lem when compar ing with the observations (e.g., 
iCavrel et ani2ci04[ ). but now is consistent if NLTE 
effects are taken into account, as shown in K06. 

We now compare these model calculations with the 
available abundance measur ements in the extr emely 
metal- poor DLA reported by ICooke et al.l H2010bl . filled 
circles). This DLA was originally identified in the 
Sloan Digital Sky Survey (SDSS) spectrum of the QSO 
J0035— 0918. Follow-up high resolution spectroscopy in- 
dicated Zabs = 2.3400972, log iV(HI) /cm"^ = 20.55±0.1, 
[Fe/H] ~ —3.04, and a pronounced carbon enhancement 
[C/Fe] ~ -KL53. 

• The observed [C/Fe] ratio is well reproduced with 
our faint SN/HN models. Low- mass AGB stars 
(1 — 4Mfyi) can also provide such high [C/Fe] (e.g., 
iKarakasI 120101 ). However, such low- mass stars are 
unlikely to contribute at the redshift of the C- 
rich DLA. With the star formation history in the 
solar neighborhood, [C/Fe] reaches the maximum 
value at z = 1.8 due to the AGB contribution 
(|Kobavashi et al.|[20Tol) . 

• The [O /Fe] ratio is higher than for non- faint super- 
novae ([0/Fe] - 0.5 - 0.6 for 25Mq in K06) and is 
consistent with the faint SN model because of the 
smaller ejected Fe mass. 

• The [(Si,S)/Fe] ratios are similar to those of non- 
faint supernovae, and also consistent with the faint 
SN/HN models. 

• The low Al abundance strongly suggests that the 
enrichment source is not Pop H supernovae but pri- 
mordial supernovae. Even for Z — 0.001, [Al/Fe] 
is larger than 0.35 for > 2OM0 (K06). 

• The low N abundance is also consistent with the 
faint SN/HN models, and the high [C/N] ratio 
cannot be ex plained with mass lo ss from rotating 
massive stars (|Mevnet et al.|[20Tol) or intermediate- 
mass AGB stars (> 4M©). 



In order to discuss the detailed explosion mechanism 
of the supernova, it is necessary to obtain the elemen- 
tal abundances of iron-peak elements. For the C-rich 
DLA, because of the low metallicity, it is impossible 
to detect heavier elements than S except for Fe. How- 
ever, the abundance profiling approach can be used for 
other absorption systems with peculiar abundance pat- 
terns, because such systems have presumably been en- 
riched by only a small number of supernovae. In Fig- 
ure [1] we overplot (open circles) the element abundances 
recently reported by Cooke et al. (2010a) for the Zabs = 
1.62650 DLA in front of the gravitationally lensed quasar 
UM637A. This DLA, which has logiV(HI)/cm-2 = 
20.7 ±0.1 and [Fe/H] = -1.56 ±0.03, exhibits a different 
abundance pattern from other DLAs with similar metal- 
licities and also from the average population of Galactic 
metal-poor stars. In particular, Ti, Ni, and Zn are defi- 
cient relative to Fe (see Fig. 11 and 12 in ICooke et al.l 
[2010a). 

We show that this abundance pattern can also be ex- 
plained with faint supernovae. The normal [(Si,S)/Fe] 
ratios are consistent with the supernova scenario. The 
low Al abundance strongly suggests that the enrichment 
source is primordial supernovae. The observed [Zn/Fe] 
ratio is more consistent with the faint SN model (solid 
line) than the faint HN model (dashed line) or non-faint 
supernovae without mixing-fallback ([Zn/Fe] —1.7 in 
K06). This is also supported by the observed [Ni/Fe]. 
[Cr/Fe] is consistent with both supernova models. We 
should note that there is a problem in the Ti nucle- 
osynthesis yields (K06), and we do not include Ti in 
the abundance profiling. At this metallicity, however, 
there is no information on C enhancement in the DLA 
observations since C line is saturated. Therefore, it is 
uncertain how efficient the mixing-fallback is. In fact, 
the observed abundances can be explained with medium 
mixing-fallback of non-faint supernovae. In this case, the 
ejected Fe mass is larger, which may be consistent with 
the observed metallicity, and thus [C/Fe] is lower than 
for the C-rich DLA. 

In Figure [1] the dotted line is for the nucleosynthesis 
yields of the PISN o f a 17 0M(7) star, which are taken from 
lUmeda fc Nomotol (120021 ). The primordial stars with ini- 
tial masses of ~ 140 — 270Mq enter into the electron- 
positron pair instability region during the central oxygen- 
burning stages and contract quasi-dynamically. Then the 
central temperature increases, central oxygen burning 
takes place explosively, and the generated nuclear energy 
is large enough to disrupt the stars completely without 
leaving compact remnants. Compared with core-collapse 
supernovae, the abundance pattern of PISNe can be sum- 
marized as follow. 

• The odd-Z effect is much larger than ^ 1 dex. 

• [(Si, S, Ar, Ca)/Fe] are much larger than 
[(0,Mg)/Fe] because of more extensive explo- 
sive oxygen burning. 

• [Cr/Fe] is much larger because of the larger incom- 
plete Si-burning region. 

• [(Co, Zn) /Fe] are much smaller because of the much 
larger ratio between the complete and incomplete 
Si-burning regions. 
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All of these characteristics disagree with the observed 
elemental abundances of these metal-poor DLAs. Even 
at high-redshift, there is no signature of the existence 
of PISNe. The [Si/C] for PISNe is as large as +1.5, 
which is also inconsistent with the observational esti- 
mate in the interga lactic medium (IGM) ([Si/C] ^ 0.77, 
lAguirre et alll2004[ ). The IGM abundance looks more 
consistent with normal (non-faint) core-collapse super- 
novae with [C/Fe] ~ and [Si/Fe] - 0.7 (K06). 

3. DISCUSSION 

We have shown that the enrichment source of the ex- 
tremely metal-poor DLA is very likely a primordial su- 
pernova that is faint as a result of mixing-fallback to 
form a 3 — GM© black hole. It is interesting that the 
observed DLA abundance is very similar to those of 
EMP stars in the solar neighborhood including the ultra 
metal-po or star HE0557-48 40 ([Fe/H] = -4.7 5, [C/Fe] 
= +1.6. iNorris et al.ll20?)7h and BD+44°493 l|Ito et al.l 
l2009f ). Without rotation, the chemical enrichment from 
very massive stars (>50 — lOOM©) is very small because 
the central part is supposed to fall onto the black hole. 
In the early stages of galaxy formation, chemical enrich- 
ment is likely to be driven by core-collapse supernovae 
from ~ 20 — 5OM0 stars, although the mass range de- 
pends on the rotation. 

How did the first supernova enrich the first galaxy? 
In our supernova scenario, we assume that the first su- 
pernova occurs in a primordial gas cloud with the to- 
tal r nass of 10^~^M©. In hyd rodynamical simulations 
('e.g.. lMac Low &: Ferraralll999fl . the interstellar medium 
(ISM) is ionized by the supernova explosion, and the 
HI mass rapidly decreases. Afterwards, the HI mass 
slightly increases due to the recombination, and then de- 
creases due to galactic mass-outflow. Some of the ISM 
could have the HI column density of logA^(HI) ^ 20.5 
and could be observed as DLAs. The observable HI re- 
gions could have Af(HI) BOOOMg, ahhough the HI 
mass highly depends on the total mass of the galaxy, 
the radial density profile, and the inhomogeneity of the 
ISM. With our faint supernova models, the ejected C 
mass is ~ O.2M0, which is roughly consistent with 
the observed C-rich DLA. In iCooke et"all (|2010b[ ). the 
masses of carbon and neutral gas of the C-rich DLA 

are estimated as M(CII) ~ 2 (n(H)/l cm'^) Mq and 
Af(HI) - 2.5 X 10^ (rt(H)/l cm^^y^ Mq, which are con- 
sistent with our scenario if n(H) ^ 3cm^^. 

Are there any signatures of the first stars in other 
metal-poor systems? For globular cluster systems 
(GCSs), the present stellar mas s and half-light radiu s are 
10'^-^ Mq and 1 - 35 pc (e.g., IGilmore et all[2007[ ). re- 
spectively, which imply high densities {n ~ 1000). Since 
the contribution of dark matter is small, the total mass 
would be lO^^^M© with the range of star formation ef- 
ficiencies (0.001 — 0.1). The progenitors of GCSs could 
be more massive if the GCS have lost a large fraction of 
stars by the relaxat ion as shown in N-body simulations 
(jLamers et al.ll2010l ). Star formation takes place quickly, 
and > 10 — 1000 supernovae are expected to occur from 
the stellar mass with the Salpeter initial mass function, 
which might blow away the ISM and quench star for- 



mation. This is consistent with the narrow metallicity 
distribution functions and the lack of the scatt er in el- 
emental abundance ratios ([Carretta et al.|[2009[ ). In the 
elemental abundance patterns, no carbon enhancement 
that seems to be originated from supernovae is seen. 

For dwarf spheroidal galaxies (dSphs), the present 
stellar mass and half-light radius a re 10^ "^M© and 
20 - 1000 pc fe.g.. IGilmore et al.l[2007l ). respectively, im- 
plying much lower densities than in GCSs. The dark 
matter mass is about lO^M© independent of the stellar 
mass (Geha et al. 200_9). From the stellar mass, 1 — 10** 
supernovae are expec ted. As seen in the observed color- 
magnitude diagrams (|Tolstov et al.ll2009[ ). star formation 
takes place slowly with a very low rate, and thus the 
ISM is likely to be inhomogeneous. The elemental abun- 
dance pattern can be used for the abundance profiling. 
In fact, a few EMP stars at [Fe/H] < - 3.5 in dSphs 

do show carbon enhancement ()Norris et al.ll2010[ ). which 
suggests a contribution from faint supernovae. Also in 
the outer halo of the Milky Way Galaxy, three stars at 
[Fe/H] < - 2.5 do show [C/Fe]' - +2 (Beers 2010, pri- 
vate communication); such stars may be disrupted from 
dSphs. 

4. CONCLUSIONS 

In the early stages of chemical enrichment, the inter- 
stellar medium is supposed to be highly inhomogeneous, 
so that the properties of the first objects can be directly 
extracted from the comparison between the observed ele- 
mental abundances and nucleosynthesis yields. We have 
shown that the observed abundance pattern of the very 
metal-poor C-rich DLA is in excellent agreement with the 
nucleosynthesis yields of a primordial star that explodes 
as a faint core-collapse supernova owing to the efficient 
mixing and fallback. The nucleosynthesis yields of PISNe 
are not consistent with the observation. The contribution 
from rotating massive stars seems to be small because of 
the lack of N enhancement. The contribution from AGB 
stars should be very small because of the N abundance 
and of the enrichment timescale. Since the DLA abun- 
dances reflect the chemical enrichment in gas-phase, the 
binary or accretion scenarios of the EMP stars do not 
work. Thus, we conclude that enrichment by primordial 
supernovae is the best solution to explain the abundance 
pattern of the C-rich DLA. The abundance pattern of the 
C-rich DLA is similar to those of EMP stars such as the 
ultra metal-poor star HE0557-4840. Some of EMP stars 
in dSphs and the Galactic outer halo also show similar 
carbon enhancement at [Fe/H] < — 3. Chemical enrich- 
ment by the first stars in the first galaxies is likely to be 
driven by core-collapse supernovae. 
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